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Abstract 

The phenomenon of frequency and phase synchronization in stochastic systems requires a revision 
of concepts originally phrased in the context of purely deterministic systems. Various definitions 
of an instantaneous phase are presented and compared with each other with special attention 
payed to their robustness with respect to noise. We review the results of an analytic approach 
describing noise-induced phase synchronization in a thermal two-state system. In this context 
exact expressions for the mean frequency and the phase diffusivity are obtained that together 
determine the average length of locking episodes. A recently proposed method to quantify frequency 
synchronization in noisy potential systems is presented and exemplified by applying it to the 
periodically driven noisy harmonic oscillator. Since this method is based on a threshold crossing 
rate pioneered by S.O. Rice the related phase velocity is termed Rice frequency. Finally, we discuss 
the relation between the phenomenon of stochastic resonance and noise-enhanced phase coherence 
by applying the developed concepts to the periodically driven bistable Kramers oscillator. 
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Studying synchronization phenomena in stochastic systems necessitates a re- 
vision of concepts originally developed for deterministic dynamics. This state- 
ment becomes obvious when considering the famous phase-locking effect pi |2[: 
unbounded fluctuations that occur, for instance, in Gaussian noise will always 
prevent the existence of a strict bound for the asymptotic phase difference of 
two systems. Nevertheless, a reformulation of the synchronization phenomenon 
in the presence of noise is possible by quantifying the average duration (Ti oc k) of 
locking epochs that are disrupted by phase slips. In case that (Xi oc k) ^> T , where 
To is some characteristic time of the dynamics, e.g. the period of an external 
drive or the inverse of some intrinsic natural frequency, it is justified to speak 
about effective phase synchronization. 



I. INTRODUCTION 



From the conceptual point of view different degrees of synchronization can be distin- 
guished: complete synchronization ||, generalized synchronization ||, lag synchronization 
f|, phase synchronization || [7j, and burst (or train) synchronization ||. In the following, 
we focus our attention on phase synchronization in stochastic systems that has attracted 
recent interest for the following reason: in many practical applications the dynamics of a 
system, though not perfectly periodic, can still be understood as the manifestation of a 



stochastically modulated limit cycle [§, |10| . As examples, we mention neuronal activity |]TT 
the cardiorespiratory system ||12|| , or population dynamics fl3| . 

Given a data set or some model dynamics there exists a variety of methods to define an 
instantaneous phase <j)(t) of a signal or a dynamics. For a clear cut separation of deterministic 
and noise-induced effects it is essential to assess the robustness of each of these different 
phase definitions with respect to noise. Section II is devoted to this issue. Since we do not 



distinguish between dynamical and measurement noise our treatment is also tied to the 
question how synchronization can be detected within any realistic experimental data. 

The synchronization properties of a noisy system can be classified in a hierarchical man- 
ner: stochastic phase locking always implies frequency locking while the converse is not 
true in general. On the other hand, small phase diffusivity is necessary but not sufficient 
for phase synchronization. This will become clear in Sec. fill when we review an analytic 
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approach |TJj to stochastic phase synchronization developed for a thermal two-state system 
with transitions described by noise-controlled rates. 



A recently proposed method |I6| to measure the average phase velocity or frequency in 
stochastic oscillatory systems based on Rice's rate formula for threshold crossings [IT], |TB[ will 
be presented and discussed in Sec. [TV|. The Rice frequency proves to be useful especially in 
underdamped situations whereas the overdamped limit yields only finite values for coloured 
noise. Its relation to the frequency based on the widely used Hilbert phase (cf. Sec. [II C|) is 
discussed and illustrated. 

In the final Sec. [V| we connect the topic of stochastic resonance (SR) [pi], |2(| with results 
on noise-enhanced phase coherence [21]. To this end we study the synchronization properties 



of the bistable Kramers oscillator driven externally by a periodic signal p2| . As a comple- 



ment to the frequently investigated overdamped limit, we consider here the underdamped 
case employing the methods presented in Sec. [TV|. 

II. PHASE DEFINITIONS IN THE PRESENCE OF NOISE 
A. Natural phase <p N 

A phase occurs in a quite natural way when describing the cyclic motion of an oscillator 
in phase space. Self-sustained oscillators are nonlinear systems that asymptotically move on 
a limit cycle. The instantaneous position in phase space can be represented through instan- 
taneous amplitude a N (t) and phase <fi N (t). A systematic approach to relate the amplitude 
and phase dynamics to the dynamics formulated in original phase space was developed by 
Bogoliubov and Mitropolski ||23|| . Their method starts from the following decomposition of 
the dynamics 

x = v , (1) 



v 



w x + f(x,v,t,£,...) (2) 



where the function / comprises all terms of higher than first order in x (nonlinearities), 
velocity dependent terms (friction), and noise. In their work Bogoliubov and Mitropolski 
considered the function / to be a small perturbation of order e which means that the system 
is weakly nonlinear and the noise or the external forces are comparatively small as not to 



3 



distort the harmonic signal too much. The definition of an instantaneous phase proceeds by 
expressing the position x and the velocity v in polar coordinates a N and 4> N 



x (t) = a N (t) cos [(f> N (tj\ , 
v(t) = -too a N (t) sin [<f> N (t)] 



which yields by inversion [24 



it) = sfxHt) + [v(t)/u ] 

V{t)/UJQ 



arctan 



(3) 
(4) 



(5) 
(6) 



x(t) J 

It should be noted that a meaningful clockwise rotation in the x, f-plane determines angles 
to be measured in a specific way depending on the sign of u Q . Using Eqs. @, (|), @ and 
(§) it is straightforward to transform the dynamics in x and v, Eqs. ([!]) and (||]), into the 
following dynamics for a N and <p N 0, |25| 

/ (a N cos( 



a N 



[a cos 



-u a N sm( ( j) N ),t^) 



iv 

), — ujq a" sm{ 



sin 



/ JV 



N 



COS 



)■ 



(7) 
(8) 



The line x = corresponds to angles <p N = it/ 2 + mr, n G AT. As can be read off from 
Eq. (H), the phase velocity always assumes a specific value for x = fl26j , i.e., 



(x = 0) = uj . 



(9) 



This has the following remarkable consequence. We see that even in the presence of noise 
passages through zero in the upper half plane v > are only possible from x < to x > 0, in 
the lower half plane only from x > to x < 0. This insight becomes even more obvious from 
a geometrical interpretation: as the noise exclusively acts on the velocity v, cf. Eq. (0), it 
can only effect changes in the vertical direction (in x, v -space). Along the vertical line x = 0, 
however, the angular motion possesses no vertical component while radial motion is solely 
in the vertical direction and, therefore, only affected by the noise. From this we conclude 
that between subsequent zero crossings of the coordinate with positive velocity the phase 
has increased by an amount of 2n. This finding establishes a simple operational instruction 
how to measure the average phase velocity of stochastic systems. We will come back to this 



point in Sec. IV 



B. Linear interpolating phase (f> L 



As we have just seen zero crossings can be utilized to mark the completion of a cycle. 
This can be generalized to the crossings of an arbitrary threshold with positive velocity or 
even to the crossing of some separatrix. In this connection the concept of isochrones of a 
limit cycle has to be mentioned [E7|. All of these extensions of the natural phase require 



a thorough knowledge of the dynamics and the phase space structure. In many practical 
applications, however, the detailed phase portrait is not known. Instead, one is given a data 
series exhibiting a repetition of characteristic marker events, e.g. the spiky peaks of neural 
activity, the R-peaks of an electrocardiogram, or pronounced maxima as found in population 
dynamics. These marker events can be used to pinpoint the completion of a cycle, k, and 
the beginning of a subsequent one, k + 1. It is then possible to define an instantaneous phase 
4> L (t) by linear interpolation, i.e., 

<P L (t) = + t ~ t \ 2tt + k2n (t k <t< t k+1 ) (10) 

where the times t k are fixed by the marker events. Reexpressing the time series x(t) of the 
system as 

x(t) = a L (t) cos[0 L (t)], (11) 

then defines an instantaneous amplitude a L (t). The benefit of such a treatment is to reveal 
a synchronization of two or more such signals: whereas the instantaneous amplitudes and, 
therefore, the time series might look rather different, the phase evolution can display quite 
some similarity. If the average growth rates of phases match (notwithstanding the fact that 
phases may diffuse rapidly) the result is termed frequency locking. Small phase diffusion, in 
addition to frequency locking, means that phases are practically locked during long episodes 
that occasionally are disrupted by phase slips caused by sufficiently large fluctuations. This 
elucidates the meaning of effective phase synchronization in stochastic systems. 

As should be clear from its definition the linear phase relies on the clear identification of 
marker events. With increasing noise intensity this identification will fail since sufficiently 
large fluctuations may either mask true or imitate spurious marker events. On the other 
hand, in some cases, e.g. for excitable systems, fluctuations can be essential for the generation 
of marker events, i.e., marker events may be noise- induced. 
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As a final remark, let us mention that relative maxima of a differentiable signal x(t) 
correspond to positive-going zeros of its derivative x(t). However, this seemingly trivial 
connection is overshadowed by complications if the derivative itself is a non-smooth function 
that does not allow to easily extract the number of zero crossings (cf. Sec. [TV] and Fig. |5|). 



C. Hilbert phase 4> H 

In situations where a measured signal x(t) exhibits a lot of irregularity it is not quite clear 
how to define a phase - the signal might look far from a perturbed harmonic or even periodic 
one and marker events cannot be identified unambiguously. The concept of the analytic 
signal as introduced by Gabor |Q offers a way to relate the signal x(t) to an instantaneous 
amplitude a H (t) and a phase <fi H {t). The physical relevance of a such constructed phase 
is a question of its own; for narrow-band signals or harmonic noise it has a clear physical 
meaning whereas the general case requires further considerations (cf. Appendix A2 in ||29|j). 

The analytic signal approach extends the real signal x(t) to a complex one z(t) = 
x(t) + iy(t) = a H (t) exp[(p H (t)} with the imaginary part y(t) resulting from an appropri- 
ate transformation of the real part x(t). Instead of taking y(t) = —x(t)/cuo as for the 
natural phase we search y(t) as the result of a convolution of x(t) with some appropriate 
kernel K(t), i.e., 

oo 

y(t)= J x{r)K{t - t) dr . (12) 

— oo 

Now, appropriate means that the kernel has to be chosen such that the method reproduces 
the phase of a harmonic signal. Applying the convolution theorem it is easy to see that 
the Fourier transform Y(uj) of the transformed signal y(t) should be related to the Fourier 
transform X(u) of the original signal x(t) by a phase shift that transforms a cosine into a 
sine, i.e., 

Y(u) = -isgn(co)X(uj) . (13) 

where sgn(-) is the sign function. By an inverse Fourier transform we thus find that K(t) = 
l/(7rt) which implies that y(t) is related to x(t) via the Hilbert transform 

oo 

y(t)=x H (t) = -P [ p^alr. (14) 

7T J t — T 
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The symbol P in front of the integral in Eq. ([TJ|) is a reminder that the integral has to be 



evaluated in the sense of the Cauchy principal value. The fact that the Hilbert phase |p4 

~x H (t) 



4> H {t) = arctan 



(15) 



_ x(t) 

arises as the result of a convolution instead of a differentiation makes it less sensitive to short- 
lived small fluctuations. This observation was already reported by Vainstein and Vakman 

5g. 

Moreover, the construction by a convolution brings the Hilbert phase in close contact 
with the wavelet transform that is widely used to compute a time-dependent spectral de- 
composition of non- stationary signals ||31|| . By virtue of Eq. (|l"3|) it is evident that the Fourier 
transform of the analytic signal Z(u) satisfies the relation 

Z(u) = [1 + sgn(cj)] X{u) = 2 9(u) X{u) (16) 

where O(-) is the step function. Equation ([IB]) shows that all positive Fourier components 
of the original signal contribute with equal weights. Selecting, instead, a certain frequency 
band by using a subsequent Gaussian filter (in frequency space) corresponds to a convolution 
with the Morlet wavelet (Gabor function). Generalized phase definitions and the wavelet 
transform were employed to detect phase-synchronous activity in the brain and in a 
chaotic laser array 



D. Discrete phase 



lD 



Multistability is one of the crucial consequences of nonlinearity and plays a dominant 
role for many important topics, e.g. evolution, information processing and communication, 
pattern formation, etc. Frequently, fluctuations play a benificial role in that they effect 
transitions between the different states that are directly tied to the performance of a task. 
The phenomenon of SR |I9|, |20||, for instance, can be observed in a bistable system. Typically, 
the information that is processed in a bistable system does not require to keep track of a 
continuum but is rather contained in the switching events between the two states |34| . Hence, 



it is desirable to link the dichotomous switching process to a description in terms of an 
instantaneous phase. Since switching to and fro constitutes one cycle and, thus, corresponds 
to a phase increment of 2tt the linear interpolating phase <fi L can be readily constructed 
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FIG. 1: The linear interpolating phase, the Hilbert phase and the discrete phase agree in the 
description of a dichotomous switching process 

employing switchings as marker events of half-cycles. Furthermore, the Hilbert phase can 
be easily computed. Alternatively, it is possible to use the switching process (between — 1 
and +1) and construct a discrete instantaneous phase <j) D {t) changing discontinuously at the 
switching events. The last-mentioned instantaneous phase is obtained simply by multiplying 
the number of switches (or the number of renewals in renewal theory ||35|| ) k(t) by the value 
7r, i.e., 4> D {t) = k(t)iT. The instantaneous state, in turn, can be obtained from the discrete 
phase via x(t) = exp[i<f) D (t)]. In Fig. [I] we show how the three alternative phases <fr L ,<f) H 
and <p D agree in the description of a dichotomous switching process. Note that the natural 
phase 4> N it) is related to the underlying process in real phase space and, hence, cannot be 
deduced from the two-state signal. The advantage of the discrete phase is that it allows an 
analytic treatment of effective phase synchronization in stochastic bistable systems. This 



will be addressed in Sec. Ill 



The robustness of the discrete phase with respect to noise comes into play not when mak- 
ing the "transition" from the switching process to the instantaneous phase <j> D (t) but when 
constructing the switching events from the continuous stochastic trajectory. At the level 
of a Markovian switching process noise enters only via its intensity that changes transition 
rates. 



S 



E. Phases in higher dimensional systems 

In systems of dimension d > 1 there are various ways how to define one or even more in- 
stantaneous phases: projecting the 2c?-dimensional phase space onto 2-dimensional surfaces, 
choosing Poincare sections or computing the Hilbert phase for each of the coordinates. Many 
of the methods can only be done numerically and always require to consider the dynamics in 
detail in order to check whether a made choice is appropriate. We will not elaborate these 
details here but refer to J^J (especially Chap. 10) and references therein. 



III. ANALYTIC TREATMENT OF A DRIVEN NOISY BISTABLE SYSTEM 



A. Setup of the doubly dichotomous system 

In this section we consider a stochastic bistable system, for instance a noisy Schmitt 
trigger, which is driven externally either by a dichotomous periodic process (DPP) or a 
dichotomous Markovian process (DMP). The dichotomous character of the input shall be 
either due to a two-state filtering or be rooted in the generation mechanism of the signal. 
The bistable system generates a dichotomous output signal. For convenience we choose to 
label input and output states with values +1 and —1 respectively. The DPP is completely 
specified by its angular frequency Q = vr/T where T is the half-period. Accordingly, the 
DMP is fully characterized by its average switching rate 7 = Tq . 

Since both input and output are two-state variables it is possible to study phase synchro- 



nization in terms of discrete |5T| input and output phases <j>{ n (t) and <fi ut(t) respectively. 
Consequently, also the phase difference (f>(t) = 4> out (t) — ^> in (t) is a discrete quantity that 
can assume positive and negative multiples of it. From the definition of the phase difference 
<p it follows that each transition between the output states increases (f) by n whereas each 
transition between the input states reduces <fi by 7r. 

Transitions between the input states are governed by the rates 

00 

W2 MP = 1 and W° pp (t)= S(t-t n ) (17) 

n=— oo 

where a single realization of the DPP is characterized by deterministic switching times 
t n = (mr + ip )Q~ 1 . Here, <^ is the initial phase of the input signal rendering the periodic 
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W™\<t>) = g{(f>) = a exp 



(18) 



process non- stationary (cyclo-stationary). To achieve strict stationarity we average the 
periodic dynamics with respect to y?o which is equidistributed over the interval [0,2-71"). 

In the absence of an input signal the two states are supposed to be symmetric and the 
hopping rates for both directions are identical and completely determined by a prefactor ao, 
the energy barrier AU, and the noise intensity D. The central assumption of our analysis 
is that the input signal modifies the transition rates of the output solely through the phase 
difference <fi in the following way 

AU + Aa((j)) ~ 
D 

where the function <r(0) = cos(0) = ±1 and the amplitude A < AU = 0.25 to keep the 
signal subthreshold. This definition introduces two noise-dependent time scales 

ai = a(D) exp ^— — ^ and a 2 = a(D) exp (^j^) (19) 

with a(D) = aoexp(—AU/D). The function <j(4>) favours phase differences with even 
multiples of n, i.e., in-phase configurations. 

A description of the stochastic evolution of the phase difference is based on the prob- 
abilities P(4>,t\<j)Q,to) to experience a phase difference at time t conditioned by a phase 
difference 0o & t time to- Due to the discrete character of <fi (allowing only for multiples of 7r) 
we briefly denote Pk = P{4> = kn,t\<j)o,to). Then the probabilistic evolution operator reads 
with <7k = g((f) = kir) from (|T8| ) 

dP (t) 

= LP k (t) + g k ^P k ^(t) - g k P k (t) . (20) 

While the last two terms on the right hand side of Eq. (p0|) account for the change of (f) due 
to transitions of the output the operator L reflects switches of the input 

LP k = W ± (P k+1 - P k ) (21) 

with the related input switching rates W± given by Eq. (0). As mentioned above the non- 
stationary (cyclo-stationary) character of the DPP can be cured by averaging over the initial 
phase ip . Since "temporal" and "spatial" contributions in Eq. (|21| ) are separable we can 
perform this average prior to the calculation of any moment of <fi yielding 

(^"W^E^-^H- (22) 

n n=— oo v ' 
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From Eq. fl22|) we see that the ^-averaged DPP formally looks equivalent to a DMP with 
the transition rate VL/tt. Of course, initial phase averaging does not really turn a DPP into 
a DMP. The subtle difference is that while 0-moments of the DMP continuously change in 
time related ^-averages of the DPP (before the (^-average) are still discontinuous, hence, 
temporal derivatives of functions of (0) have to be computed with care before initial phase 
averaging [IIR KHj. 



B. Noise-induced frequency locking 



Using standard techniques |37j we can derive the evolution equation for the mean phase 
difference (0) from Eq. (|20|) 



(0) = ~{u- m ) + (uj out ) 



7T 



-(^ in ) + ^ (»i + a 2 ) 



7T 



(a 2 - a 1 )(a) 



(23) 
(24) 



Here, (ui n ) denotes the average frequency of the input phase and equals 'jtt for the DMP 
and Q for the DPP. Assuming higher moments uncoupled, i.e., (cx(0)) oc <r((0)), Eq. (|23| 



is Adler's equation [38| arising in the context of phase locking. Note that here both the 
frequency mismatch A = — (ui n ) + ^(ai + a 2 ) and the synchronization bandwidth A s = 
^(a 2 — ai) are noise dependent. This elucidates the opportunity to achieve noise-induced 
frequency and effective phase locking. For the short-time evolution a necessary condition 
for locking is |A| < A s which defines "Arnold tongues" ]l| of synchronization in the A vs. D 
plane. 

The kinetic equation for (a) can be evaluated explicitly yielding 



2- h ai + a 2 



7T 



(a) + a 2 — ai . 



From Eq. 



we see that (a) approaches a stationary value 

a 2 — ai 



2 h ai + a 2 



(25) 



(26) 



7T 



that exactly coincides with the stationary correlation coefficient between the input and 
output [f4l] . The relaxation time is given by r = [2(u;i n )/7r-|-ai-|-a 2 ] _1 . Hence, the stationary 
output phase velocity can be achieved from Eq. ( ]23|) by insertion of Eq. (|26"D yielding 

(^out) = + a 2) - ^ ( fl 2 - ai)<<7*) . (27) 
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FIG. 2: Mean output switching frequency (in units of n) vs. noise intensity for the DMP [or the 
DPP] signal for 7 = 0.001 [or £l/ir = 0.001 respectively] and for three amplitudes A: (solid), 0.1 
(dotted), 0.2 (dashed) (AU = 0.25, other parameters see text). "Arnold tongues", defined by the 
demand of "sufficiently small slopes" (see text), are shown in the inset for various 7: 0.001 (solid), 
0.002 (dotted), 0.005 (dashed). 



This expression is in agreement with similar results derived in the context of resonant acti- 
vation p9| . In Fig. H we depict the mean output switching rate (ui* ut ) as a function of the 
noise intensity D for several values of the input signal amplitude A. With increasing ampli- 
tude the region of frequency locking, i.e., the region of D for which (ui* ut ) ~ (u>i n ), widens 
(cf. also Fig. 9 in [fHJ). For most of these intensities the bistable system possesses rates 
that do not obey the time scale matching condition; nevertheless, on the average the output 
switching events are entrained by the input signal. As mentioned before the whole effect is 
non-linear: If detuning becomes too large the output gets desynchronized and returns to its 
own dynamics. 

Besides many alternatives, we denned the width of the frequency locking region by requir- 
ing the slope of the curves to be less than 30% of the slope for A = at the point where the 
output switching rate coincides with 7 [or Q/n respectively] - simultaneously disqualifying 
the initial flat region for small D. The resulting "Arnold tongues" (compatible with data 
from [21, 41, 42]) are shown in the inset of Fig. ||. It can be seen that frequency locking 



necessitates to exceed a minimal amplitude that shifts to lower values for slower signals [42 



Let us emphasize that the frequency locking region groups around the noise intensity that 
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satisfies the time scale matching condition (in our case by definition u in = w* ut (£))); for a 
harmonic input this range of D also maximizes the spectral power amplification, in contrast 
to values of D where the signal-to-noise ratio attains its maximum (D oc AU) [19|, [2(J . 



C. Phase locking and effective synchronization 

The phenomenon of phase locking can be demonstrated by considering the diffusion 
coefficient T> of the phase difference, achieved as the time derivative of the variance |<9 t [(0 2 ) — 
(0) 2 ]. Performing the calculation for both the DMP and DPP yields 



7T 



v = v in + — 



(^out) 



7T 



a 2 - a lj 



71 



(28) 



with D in = ^-7 for the DMP and V in = for the DPP. The stationary correlator ((8(f)8a)*), 
i.e., the asymptotic limit of (4>a) — ((f)) (a), can be computed from the corresponding kinetic 
equation ]T5[. Inserting ((5<fiSa)*) into Eq. Q2lf ) we thus find for the DMP (cf. Fig. || top) 



V 



DMP 



71 

y 



7+ 



*\2 



71 



- (2 7 - ( ai + a 2 )) (a*) 



(29) 



and for the DPP (cf. Fig. | bottom) 



V 



DPP 



TT 

T 



71 



2--(ax + a 2 ) ) (a*) 

TT 



*\2 



-^(«2-a 1 )(a*)(l + (a*) 2 ) + -( ( x*) 

Z 71 



(30) 



with (a*) given by Eq. @ and (a;* ut ) by Eq. (0). Both Eqs. (f2p]) and (§Dp possess the 
same structure V = V in + V out — V co with V out = ~ (o;* ut ). Since is never decreasing 
(with increasing noise intensity) the same is true for the sum of the first two terms. The 
possibility of synchronized input-output jumps is rooted in T> co . Since this term comprises 
only contributions scaling with powers of (ex*), which itself rapidly vanishes for small D 
(cf. inset of Fig. |] bottom), we first observe an increase of T>. An increase of (a*) signals the 
coherent behaviour of input and output and, consequently, endows V co with considerable 
weight to outbalance the increase of V out . As can be seen from the inset of Fig. |] bottom 
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FIG. 3: Effective diffusion coefficients D (top) and T> (bottom) of the instantaneous phase 
difference (j> as functions of the noise intensity D for 7 = — = 0.001 and for three amplitudes A: 
(solid), 0.1 (dotted), 0.2 (dashed) (other parameters as in Fig. §). The values at D = are 
determined solely by the input diffusion and, hence, vanish for strictly periodic signals. Defining 
the region of phase locking by the demand that D DMP < T>® MP yields the tongues depicted in 
the upper inset which, as in Fig. [2], reveals a critical amplitude varying with 7: 0.001 (solid), 0.002 
(dotted), 0.005 (dashed). As can be seen from the lower inset phase locking occurs for considerably 
large (a*) only (curves for amplitudes A: (solid), 0.1 (dotted), 0.2 (dashed)). 

a negative slope is initiated only for a sufficiently large (er*). However, the range of high 
input-output correlation (a*) does not determine the range of low diffusion coefficients since 
at rather high noise intensities D the output switches with a large variance and thus, finally 
dominates over the ordering effect of T> co . 
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Plotting the boundaries of the region where X> DMP < T>P defines the tongues depicted 
in the inset of Fig. [5] top. As for the "Arnold tongues" in Fig. |^ a minimal amplitude 
varies with the mean input switching rate (u)\ n ) and shifts to lower values when considering 
slower signals. It is worth mentioning that the addition of an independent dichotomous 
noise that modulates the barrier AU can drastically reduce this minimal amplitude if this 
second dichotomous noise switches faster than the external signal H3[ . 



The minimum of T> observed in the region of frequency locking can be equivalently ex- 
pressed as a pronounced maximum of the average duration of locking episodes (T lock ) . To 
show this we note that a locking episode is ended by a phase slip whenever the phase differ- 
ence has changed, i.e., increased or decreased, by the order of ir, or 

(4> 2 ) = (0) 2 (T lock ) 2 + 2D(T lock ) = vr 2 . (31) 



This quadratic equation can be solved for (Ti ock ) []36| and by inserting the noise-dependent 
expressions for = (<fi) and T> we can compute (Ti ock )/T as a function of noise intensity D 
where To is either I/7 for the DMP or tt/Q for the DPP respectively. The results for both the 
DMP and the DPP are plotted in Fig. A pronounced maximum for intermediate values of 
noise intensity clearly proves that noise-induced frequency synchronization is accompanied 
by noise-induced phase synchronization. 



IV. OSCILLATORY SYSTEMS AND THE RICE FREQUENCY 
A. General relations for potential systems 

As mentioned in Sec. II A positive-going zero crossings can be used to count completions 
of a cycle in oscillatory systems. In this view the average frequency, i.e., the average phase 
velocity, turns out to be the average rate of zero crossings which is captured by a formula 
put forward by Rice |17|, |18| . This elementary observation yields a novel way to quantify 
the average frequency of a phase evolution, henceforth termed the "Rice frequency" , and to 
prove frequency locking in stochastic systems. 

To detail our derivation of the Rice frequency in this section, we start from the following 
one-dimensional potential system 

x + jx + U'(x) = y/j^ + Fco8(Qt) (32) 
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FIG. 4: The normalized average duration of locking episodes (Ti oc k) (cf. Eq. (|3l"l)) for the DMP 
(top) and the DPP (bottom) exhibits an enormous maximum for intermediate values of noise 
intensity thus proving noise-induced phase synchronization 

subjected to Gaussian white noise £ of intensity D, i.e., 

<£(*)> = 0, m Z(s)) = 2D S(t - s) , (33) 

and being driven by the external harmonic force F cos (Of). In Fig. [5| we show a sample 
path for the harmonic oscillator 

x + 7 x + uj 2 x = yfi £ + Fcos(Ot). (34) 

where we used the friction coefficient 7 = 1, the natural frequency u = 1, and a vanishing 
amplitude F = of the external drive. As can be read off from Fig. [5], the velocity v = x 
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FIG. 5: Position x and velocity v of the undriven noisy harmonic oscillator Eq. (34) with friction 
coefficient 7 = 1, and natural frequency ujq = 1. Whereas the position x is smooth the velocity v 
is continuous but nowhere differ entiable. Counting of zero crossings is, consequently, only possible 
for the ^-coordinate. 

basically undergoes a Brownian motion and, therefore, constitutes a rather jerky continuous, 
but generally not differentiable signal. In particular, near a zero crossing of v there are many 
other zero crossings. In contrast to that, the coordinate x is a much smoother signal since 
it is determined by an integral over a continuous function 

t 

x{t) = x(0) + / v(r)dr, (35) 



and, therefore, differentiable. In particular, near a zero crossing of x there are no other zero 
crossings. In the following, we will take advantage of this remarkable smoothness property 



of x that is an intrinsic property of the full oscillatory system (|32|) and disappears when we 
perform the overdamped limit. 

In 1944, Rice 0] deduced a formula for the average number of zero crossings of a smooth 
signal like x in the oscillator equation (|32|) . In this rate formula enters the probability density 
P(x, v; t) of x and its time derivative, v — x, at a given instant t. The Rice rate for passages 
through zero with positive slope (velocity) is determined by fl8 



(f)(t) = / vP(x = 0,v;t)dv. 



(36) 
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This time- dependent rate is to be understood as an ensemble average. If the dynamical 
system is ergodic and mixing the asymptotic stationary rate (f s ) can likewise be achieved 
by the temporal average of a single realization. Let iV([0, t]) be the number of positive-going 
zeros of the signal x in the time interval [0,t]. Using ergodicity, the relation 



(fs) 



vPJx = 0, v) dv 



lim 

t— >co 



N([0,t)) 



(37) 



is fulfilled for the process characterized by the stationary density P s (x,v). In the following 
we always consider stationary quantities. As explained in Sec. II A, the zero crossings can 
be used as marker events to define an instantaneous phase <fi L (t) by linear interpolation, 
cf. Eq. (0)- The related average phase velocity is the product of the (stationary) Rice rate 
and 27r and, hence, called the (stationary) Rice frequency 



oc 

2n(f s ) = 2n J vP s (x = 0, v) dv . 



(38) 



For a dynamics described by a potential U(x) in the absence of an external driving, i.e., 
2|) with F — 0, the stationary density can be calculated explicitly yielding 

,2 



P s (x, v) = C exp 



- + U{x) ) /D 



(39) 



where C is the normalization constant. From this and the application of Eq. (|38|) , it is 
straightforward to derive the exact result 



\Jl-nD exp 



!7(0) 



D 



J exp 



U(x) 
D 



(40) 



fix- 



Without loss of generality we can set U(0) = 0. In the limit D — > 0, we can perform a 
saddlepoint approximation around the deepest minima x m (e.g. for symmetric potentials). 
In this way we find the following expression valid for D AU = U(0) — U(x,j) , i.e., the 
small noise approximation, 

r u(xa] i _1 

exp ^ L 



E 



(41) 



In the limit -D — > oo, we have to consider the asymptotic behaviour of the potential, 
\im x ^± 00 U(x), to estimate the integral in Eq. (PEJ). For potentials that can be expanded 
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in a Taylor series about zero and that, therefore, result in a power series of order 2m, i.e., 
U(\x\ — > oo) ~ x 2m , we can rescale the integration variable by x = D^l^x. For sufficiently 
large D, the integral is dominated by the power 2m term. In this way we find the large 
noise scaling 

(uj) R ~ ' D a , with a = . 42 

2m 



Applying Eqs. (|40f) and (|4lD to the harmonic oscillator (|34|) we immediately find that (o>).r = 
u>o, independent of 7 and for all values of D > 0. This is also in agreement with Eq. (ft2"l) . 
It follows because m = 1 implies that, for large noise, the Rice frequency (w)_r does not 
depend on D at all. Note, however, that in the deterministic limit, i.e., for D = 0, we have 
the standard result 

(D=o) j v^o - 7 2 /4 for 7 < 2w , , 

= \ for 7 > 2^ ' (43) 

which explicitly does depend on the friction strength 7 > 0. Therefore, the limit D = is 
discontinuous except in the undamped situation 7 = 0. 

The similarity of Eqs. fl40D and ([y]) with rates from transition state theory |44j will be 
addressed below when we discuss the bistable potential. 

B. The role of coloured noise 

It is well known that the Rice frequency cannot be defined for stochastic variables that 
integrate increments of the Wiener process (white noise). From Eq. (|32"D this holds true for 
the velocity v = x. This is so, because the stochastic trajectories of degrees of freedom being 
subjected to Gaussian white noise forces are continuous but are of unbounded variation and 



nowhere differentiable || 0. This fact implies that such stochastic realizations cross a 
given threshold within a fixed time interval infinitely often if only the numerical resolution 
is increased ad infinitum. This drawback, which is rooted in the mathematical peculiarities 
of idealized Gaussian white noise, can be overcome if we consider instead a noise source 
possessing a finite correlation time, i.e., coloured noise, see ref. |4^| . To this end, we consider 



here an oscillatory noisy harmonic dynamics driven by Gaussian exponentially correlated 
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noise z(t), i.e., 



x = v (44) 
v = —7 v — uJq x + ^7 (45) 
i = -- + -£, (46) 



r r 

with z(t) obeying (z(t)) = and 

= 7 exp(-t^) . (47) 

Following the same reasoning as before we find for the Rice frequency of x(t) as before 

00 00 



(u) x = J dv J dz vP s {0,v,z) (48) 
-00 

= • (49) 

Likewise, upon noting that within a time interval At, — At(— 7X — UqX + yFyz) < v < 0, or 
— Ai(— c^qX + ^/t 7 ^) + (9(Ai) 2 < d < 0, respectively, the Rice frequency of the zero crossings 
with positive slope of the process v(t) is given by 



00 00 



which is evaluated to read 



J dx J dz (y/j z — Uq x) P s {x, 0, z) , (50) 



-00 X 



(u) v = ^Jul + 1 -. (51) 

The result in (^) shows that for small noise colour r the Rice frequency for (u) x assumes 
a correction (u) x ~ uj ^1 — ^J, as r — > + . In clear contrast, the finite Rice frequency for 
the velocity process v(t) (f45|) diverges in the limit of vanishing noise colour proportional to 

C. Relation between Rice and Hilbert frequency 

To exemplify the relation between the Rice frequency (uj)r and the Hilbert frequency 



{uj)hi again we consider the damped harmonic oscillator Eq. (34) agitated by noise alone. 



In Fig. ^ we show a numerically evaluated sample path and the corresponding Hilbert phase 
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FIG. 6: Signal x(t) (top panel) and corresponding Hilbert phase {t)/2n modulo 1 (bottom 



panel) for the undriven harmonic oscillator Eq. (34) with friction strength 7 = 1, noise intensity 



D = 1, natural frequency ujq = 1, and driving amplitude F = 0. Note that although there are 
successive zero crossings of x with positive slope near t ~ 3 and t ~ 9 the Hilbert phase does not 
increase by 2ir. 

(normalized to 2n and modulo 1) using the parameters 7 = 1, D — 1,Uq — 1, F — 0. An 
important point to observe here is that around t ~ 3 and t ~ 9 the Hilbert phase 0^ does 
not increase by 2tt after two successive passages through zero with positive slope. This shall 
illustrate the difference between the Hilbert phase and the natural phase. In Subsec. [II C 



this observation was already mentioned as a consequence of the nonlocal character of the 
Hilbert transform. In particular, short and very small amplitude crossings to positive x 
are not properly taken into account by the Hilbert phase since they only result in a small 
reduction of (fi H . This leads us to conjecture that quite generally 

(uj)r > (u) H (52) 

holds. In fact, for the case of the harmonic oscillator that generates a stationary Gaussian 
process one even can prove this conjecture by deriving explicit expressions for (u>)r and 
(uj)h- As usual, let S(u) denote the spectrum of the stationary Gaussian process x. Then 
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the Rice frequency can be recast in the form of [IS 



R 



Juj 2 S{uj) du 



oo 

f S(u) duj 
l o 



1/2 



(53) 



A similar expression (additionally involving an Arrhenius-like exponential) exists when con- 
sidering not zero crossings, as in Eq. 



but crossings of an arbitrary threshold. In [IB it 
was shown that the Hilbert frequency of the same process x is given by a similar expression, 
namely 



f uj S{ui) dui 
o 

oo 

f S(u) du 

L 



(54) 



Interpreting the quantity S{uj)/ J °° S{ui)duj as a probability density P(uj), uj G (0, oo), we 
can use the property that the related variance is positive, i.e., 



uj 2 P{uj) du; > 



uiP{ui) dui 



Lo 



(55) 



Taking the square-root on both sides of the last inequality immediately proves Eq. (|52j). 
Using the spectrum of the undriven noisy harmonic oscillator 

4jD 



SV) 



(56) 



(u)q — uj 2 ) 2 + 7 2 cj 2 

and employing Eqs. (|53D and (^4[), it is easy to see that both (uj)r and (uj)h do not vary 
with D. We have already shown above that (uj)r = ujq. In contrast to this, (uj)h is a 
monotonically decreasing function of 7 that approaches uio from below in the limit 7 — > + . 



D. Periodically driven noisy harmonic oscillator 

The probability density of the periodically driven noisy harmonic oscillator can be de- 
termined analytically by taking advantage of the linearity of the problem. Introducing the 
mean values of the coordinate and the velocity, (x(t)) and (v(t)), the variables 

x = x — (x), v = v — (v) (57) 
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obey the differential equation of the undriven noisy harmonic oscillator. In the asymptotic 
limit t —>■ oo the mean values converge to the well known deterministic solution 

F 



<*(*)> 

(«(*)> 
8 



v /( Lc ,2_C72)2 +7 2 fi 2 



VK 2 - n 2 ) 2 + 7 2 ^ 2 



cos(fit — 5) 



sin(fit — 5) 



arctan 



7 fi 



cj 2 - fi 2 



(58) 
(59) 
(60) 



with the common phase lag S. Therefore, after deterministic transients have settled the 
cyclo- stationary probability density of the driven oscillator reads 



Pcs(x,v;t) = P s {x- (x(t)),v- (v(t))) 



with the Gaussian density 

P«(x,v) 



2ttD 



exp 



V 2 UnX 2 
h — 

2 2 



ID 



(61) 



(62) 



Using Eq. ( |38"D the cyclo-stationary probability density (|6lD yields an oscillating expression 
for the Rice frequency (u})n(t). The time dependence of this stochastic average can be 
removed by an initial phase average, i.e., a subsequent average over one external driving 
period 27r /Q, 

2tt/Q, 



R 



Qdt 



o 

2tt/Q oo 



v P cs (0, v ; t) (if OtZt . 



(63) 



(64) 



o o 



The resulting analytical and numerically achieved values of the Rice frequency as a function 
of the noise intensity D are shown in Fig. [7] for fixed u> = l,F = l,Q = 3 and various values 
of 7. For small noise intensities D the Rice frequency (u) r is identical to the external driving 
frequency Q, whereas for large noise intensities the external drive becomes inessential and 
the Rice frequency approaches (w)r = u> . 

Further insight into the analytic expression (|64]) is gained from performing the following 
scale transformations 



t = Qt-5 



and 



x 



x 



(65) 
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FIG. 7: Rice frequencies for the driven harmonic oscillator Eq. (j34|) with natural frequency ojq = 1, 
driving amplitude F = 1, and driving frequency O = 3 for different values of the friction strength 7. 
The numerically achieved values (symbols with error bars) match the analytical curves determined 
using Eq. (p4|). 



from which we immediately find the rescaled velocity 

/2D/Qdx 



dx 
dt 



2D v. 



dt 

Inserting these dimensionless quantities into Eq. ( |64D yields 



(u) R = uj I(A,u Q ) 

2tt—S 00 

I(A, uj ) = - v exp 



(v + Asintf 



-s 



(uq a cos ty 



dv dt 



where we have defined further dimensionless quantities 

~ Q F 
A = 



^0 



(66) 



(67) 



(68) 



(69) 
(70) 



Due to the 27r periodicity of the trigonometric functions, the integral (^) does not change 
when shifting the interval for the integration with respect to t back to [0, 2n]. Hence, I is 



24 



only a function of A and ujq. An expansion for small A yields 



(u) R = u Q 
which implies for large D/F 2 



1 + L^Ia 2 + 0(A 4 



(71) 



F 2 

{uj) r -ujq ~ — . (72) 

The opposite extreme, A —>■ oo or D/F 2 —>■ 0, can be extracted from a saddlepoint approx- 
imation around v = A and t = 3ir/2. Following this procedure, the integral (|68| ) gives the 
constant 1/uq. This directly implies (uj)r = fi. 

The crossover between these two extremes occurs when the first correction term in (ffll) 
is no longer negligible, i.e., for 



|l-^0l ,2 



2 

When solved for the crossover noise intensity D co /F 2 , this yields 



A « 1. (73) 



(74) 



F 2 4 [(cj 2 - ft 2 ) 2 + ( 7 fi) 2 ] ' 
which, for the parameters used in Fig. [5], correctly gives values between 10~ 2 and 10 _1 . 

In Fig. the parameters F, Q, and ojq and, hence, uj are identical for all curves. Solving 
A(ji, Di) = A (72, -D2) with respect to -D 2 shows that the curves become shifted horizontally 
as in the log-linear plot in Fig. [7[ Another way to explain this shift is by noting that 
dDcJd-f < 0. 

V. BISTABLE KRAMERS OSCILLATOR: NOISE-INDUCED PHASE COHER- 
ENCE AND SR 

A. Rice frequency and transition state theory 



The bistable Kramers oscillator, i.e., Eq. (|32j) with the double well potential 

4 2 

U(x) = - - \ , (75) 



is often used as a paradigm for nonlinear systems. With reference to Eq. ( |3"2"D the corre- 
sponding Langevin equation is given by 

x + 7X + a; 3 - x = + F cos (Vlt) (76) 
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FIG. 8: Rice frequencies for the undriven bistable oscillator Eq. ( |7q ) with friction strength 7 = 1. 
Numerical values with error bars match the analytically determined values (dotted line) using 
Eq. (|38|) with Eq. (77). As expected, for large values of D the Rice frequency scales like D 1//4 . The 
solid line presents the leading weak noise approximation in Eq. (p2|). 



which, in the absence of the external signal, F 
distribution 



0, generates the stationary probability 



1 1 y2 x ^ x<2 



/D 



(77) 



with the normalization constant C. Using this stationary probability density and Eq. ( |38| ) 
we can determine the Rice frequency analytically. In Fig. |8] we depict this analytic result to- 
gether with numerical simulation data including error bars. The simulation points perfectly 
match the analytically determined curve. As expected for the asymptotically dominant 
quartic term, i.e., m = 2 (cf. Sec. |IV A| , especially Eq. (fi2"|)), the Rice frequency scales as 
(u)r ~ Z} 1 / 4 for large values of D . 

Comparing the Rice frequency formula, Eq. (]38f) , with the forward jumping 
from the transition state theory jOJ, 



fi, TST 



Zo 1 



dx dv 8(v) S(x) v exp[—H(x,v)/D] 



where 



Z = J dxdv exp [— H(x, v)/D) 

x<0 



(78) 



(79) 
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and H(x,v) = (l/2)v 2 + (l/4)x 4 — (l/2)x 2 represents the corresponding Hamiltonian, one 
can see that the difference between both solely rests upon normalizing prefactors. Whereas 
the rate k^ ST is determined by the division of the integral Eq. (|78|) by the "semipartition" 
function Z , the rate (lj) r /2tx is established by dividing the same integral Eq. (ITS!) by the 
complete partition function Z$ 

Z Q = J dxdv exp [-H(x,v)/D] . (80) 

Particularly for symmetric (unbiased) potentials, i.e., V(—x) = V(x), this amounts to the 
relation Z = 2Z , hence, 

(u) R = nk^ ST . (81) 

At weak noise, Ef,/D » 1, this relation simplifies to 

(u) R ^^-exp[-E b /D), (82) 

wherein E b denotes the barrier height and ujq the angular frequency inside the well (uo = V2) ■ 
Indeed, in the small-to- mo derate regime of weak noise this estimate nicely predicts the exact 
Rice frequency (cf. Fig. |). 



B. Periodically driven bistable Kramers oscillator 

The periodically driven bistable Kramers oscillator was the first model considered to 



explain the phenomenon of SR |47j] and it still serves as one of the major paradigms of SR 
|T9| , [20|] . In its overdamped form it was used to support experimental data (from the Schmitt 
trigger) displaying the effect of stochastic frequency locking |21], (42| observed for sufficiently 
large, albeit subthreshold signal amplitudes, i.e., for F min < F < 2/v / 27- From a numerical 
simulation of the overdamped Kramers oscillator and computing the Hilbert phase it was 
also found that noise-induced frequency locking for large signal amplitudes was accompanied 
by noise-induced phase coherence, the latter implies a pronounced minimum of the effective 
phase diffusion coefficient 

d cS = \d t [mm - m)) 2 } m 
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FIG. 9: Numerically determined Rice frequencies of the periodically driven bistable Kramers os- 
cillator Eq. (|7^) computed with the friction coefficient 7 = 0.5 and the angular driving frequency 
= 0.01 and plotted as a function of the noise intensity D. Different curves correspond to various 
amplitudes of the harmonic drive F. For larger values of F wider regions appear where the Rice 
frequency is locked to the external driving frequency f2. 



occurring for optimal noise intensity. Based on a discrete model PSj , analytic expressions for 
the frequency and phase diffusion coefficient were derived that correctly reflect the conditions 
for noise-induced phase synchronization [|HJ for both periodic and aperiodic input signals. 

To link the mentioned results to the Rice frequency introduced above we next investigate 
the behaviour of the Kramers oscillator with non-vanishing inertia. We show numerical 
simulations for Eq. (|7(J) with the parameters Q = 0.01,7 = 0.5 and diverse values of F in 
Fig. For larger values of F, a region around D « 0.05 appears where the Rice frequency 
is locked to the external driving frequency Q. Since for larger values of the external driving 
F smaller values of the noise parameter D are needed to obtain the same rate for switching 
events, the entry into the locking region shifts to smaller values of D for increasing F. 

In Fig. [TO] we present numerical simulations for fixed F = 0.384, Q = 0.01 and different 
values of the damping coefficient 7. Note that the value of F is slightly smaller than the 
critical value F c = 2/y / 27 ~ 0.3849.... For smaller values of 7 wider coupling regions appear 
since it is easier for the particle to follow the external driving for smaller damping. 

To check whether frequency synchronization is accompanied by effective phase synchro- 
nization we have also computed the averaged effective phase diffusion coefficient, this time 
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FIG. 10: Numerically determined Rice frequency as a function of the noise intensity D for the 



periodically driven Kramers oscillator Eq. (76) with the angular driving frequency f2 = 0.01 and 
driving amplitude F = 0.384 for different values of the friction coefficient 7. For smaller values of 
7 wider regions of frequency locking appear. 



defined by the following asymptotic expression 

1 



D eS = lim - ([0(t) - <0(t)>] 2 ) 



(84) 



It should be clear that the instantaneous "Rice" phase (f>(t) was determined via zero crossings. 
The connection with the instantaneous diffusion coefficient defined in Q83| ) is established by 
applying the limit t — > 00 

t 

D cS = lim - / D eS (t) dt. (85) 



In Fig. |TT] we show numerical simulations of the effective phase diffusion coefficient D c q as 
function of noise intensity D. The phase diffusion coefficient displays a local minimum that 
gets more pronounced if the damping coefficient 7 is decreased. Indeed, phase synchroniza- 
tion reveals itself through this local minimum of the average phase diffusion coefficient D e g 
in the very region of the noise intensity D where we also observe frequency synchronization, 
cf. Fig. The qualitative behaviour of the diffusion coefficient agrees also with a recently 
found result related to diffusion of Brownian particles in biased periodic potentials |49| . A 
necessary condition for the occurrence of a minimum was an anharmonic potential in which 
the motion takes place. In this biased anharmonic potential the motion over one period 
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FIG. 11: Effective phase diffusion coefficient vs. noise intensity for the periodically driven bistable 
Kramers oscillator Eq. (j76|) with angular driving frequency O = 0.01, driving amplitude F = 0.384, 
which is close-to-threshold driving, and for different values of 7. For smaller values of the friction 
coefficient 7 phase diffusion is diminished. 
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FIG. 12: Effective phase diffusion coefficient vs. noise intensity for the periodically driven bistable 
Kramers oscillator Eq. (j76|) with friction coefficient 7 = 0.5, angular driving frequency Q, = 0.01, 
plotted for the undriven case F = and for driving with an amplitude F = 0.2. 

consists of a sequence of two events. Every escape over a barrier (Arrhenius-like activation) 
is followed by a time scale induced by the bias and describing the relaxation to the next 
minimum. The second step is weakly dependent on the noise intensity and the relaxation 
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time may be even larger then the escape time as a result of the anharmonicity. For such 
potentials the diffusion coefficient exhibits a minimum for optimal noise, similar to the one 



presented in Figs. O and 112 



The average duration of locking episodes (T lock ) can be computed by equating the second 
moment of the phase difference (between the driving signal and the oscillator) to 7r 2 [ 3~6| . A 



rough estimate, valid for the regions where frequency synchronization occurs, i.e., where the 
dynamics of the phase difference is dominated by diffusion, thus reads (Ti 0C k) = n 2 / D c r or, 
when expressed by the number of driving periods @ 



(ni 0C k) 



n 7T 



2D 



cff 



In this way we estimate from Figs, | L] and [12] (ni oc k) 
relevant D c g varying between 10 -4 . . . 10 -6 . 



150... 15000 for Q 



(86) 
0.01 and 



C. SR without noise-enhanced phase coherence 

In the previous examples we have shown how frequency synchronization, revealing itself 
through a plateau of the output frequency matching the harmonic input frequency, and 
reduced phase diffusivity together mark the occurrence of noise-enhanced phase coherence. 
Optimal noise intensities were found in the range where one also observes SR (in the over- 
damped system). In order to underline that under certain conditions SR exists but may 
not be accompanied by effective phase synchronization we present simulation results for the 
Rice frequency and the diffusion coefficient in Fig. [13] obtained for the bistable Kramers 
oscillator with a friction coefficient 7 = 1 and external frequency Q = 0.1. For noise in- 
tensities D m 0.15 the output frequency matches Q and nearby the overdamped Kramers 
oscillator exhibits the phenomenon of SR, i.e., one finds a maximum of the spectral power 
amplification [[IS, |2(J. In contrast, we neither can find a minimum in the diffusion coefficient 



nor a plateau around Q meaning that no phase coherence and not even frequency synchro- 
nization can be observed. The reason is that the external signal switches much too fast for 
the bistable system to follow; note that in the two-state description with Arrhenius rates 
the prefactor a (cf. Eq. |T9]) restricts the switching frequency from above. Noise-induced 
phase coherence requires a device with a faster internal dynamics, i.e., Q <C olq. 
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FIG. 13: For friction coefficient 7 = 1 and external driving frequency = 0.1 the bistable Kramers 
oscillator does neither exhibit frequency synchronization nor noise-enhanced phase coherence but 
still stochastic resonance occurs for noise intensities in the range of values (D ~ 0.15) where 

VI. CONCLUSIONS 

We underline that the noise-induced phase synchronization is a much more stringent effect 
than stochastic resonance. This statement becomes most obvious when recalling that the 
spectral power amplification attains a maximum at an optimal noise intensity for arbitrarily 
small signal amplitudes and any frequency of the external signal. In contrast, noise-induced 
phase synchronization and even frequency locking are nonlinear effects and as such require 



amplitude and frequency to obey certain bounds (see the "Arnold tongues" in Sec. p|) . We 
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expect that the functioning of important natural devices, e.g. communication and informa- 
tion processing in neural systems or sub-threshold signal detection in biological receptors, 
rely on phase synchronization rather than stochastic resonance. 
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